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1 Introduction 

The potential to study "new pliysics" at the next generation of accelerators 
will to a large extent depend on the ability to control the strong interaction 
background. Hence the present interest in making QCD predictions for hard 
processes as quantitative as possible, with an increasing understanding of the 
importance to study higher-twist effects which are suppressed by powers of 
the large momentum. In general, higher-twist effects reflect the "leakage" of 
contributions from large distances into the process of interest. The theoretical 
status of these corrections is well-established in the total e+e~ annihilation 
cross section and in deep inelastic scattering (DIS) (and in related quantities), 
where dispersion relations relate the physical observable to the operator prod- 
uct expansion (OPE) of a T-product of currents at small distances. The OPE 
does not allow us to calculate power corrections, but one learns from the OPE 
the particular suppression of higher twist effects - for the total e+e~ 

annihilation cross section versus 1 / for the DIS structure functions and 
the process-independence (universality) of the uncalculable higher twist matrix 
elements. The structure of power-suppressed effects in "geniune" Minkowskian 
quantities is much less understood. There are phenomenological and theoret- 
ical indications that in certain situations such corrections are large - of order 
l/Q - and numerically important at all energies available today. In this talk 
we summarize the results of Ref. 1 on the structure of power corrections and 
renormalons in Drell-Yan (lepton pair) production, and speculate whether this 
example teaches us a lesson of general validity. 

Drell-Yan production, apart from its phenomenological significance, is the- 
oretically interesting because it is the simplest hard process with two large, 
but disparate scales, if we consider the production of the Drell-Yan pair (or 
heavy vector boson) close to the kinematical threshold z ^ 1 where z = Q'^/s, 

being the mass of the pair and s the total cms energy of the colliding 
partons. In this situation gluon emission into the final state is suppressed 
by small phase space (1 — z)Q ^ Q, and causes large perturbative cor- 
rections enhanced by Sudakov-type logarithms ln(l — z). Taking moments 
a{N,Q'^) = f dzz^~^da/dz and subtracting collinear divergences by forming 
the ratio of the Drell-Yan cross section and the quark distribution squared, one 
is left with the perturbative expansion for the logarithm of the "hard" cross 
section 

Iiiludy{N, Q2) = ^^as In N + uj^al In' N + ... ujka\ \n+^ N+..., (1) 

where as = as{Q) and we have shown terms with the leading power of the large 
logarithm at each order of as- In leading logarithmic approximation (LLA), 
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these terms are resummed to all orders by the elegant formula 

\nuDy{N,Q^)^'-^ dz'- ^a.(fci). (2) 

Expansion of the running coupling under the integral in powers of as{Q^) gen- 
erates a perturbative series which correctly reproduces all coefficients in Eq. ^ 
Eq. 1^ takes into account soft and coUinear gluon emission, and can be improved 
systematically by including next-to-leading logarithms etc. The leading soft- 
coUinear contribution has a geometrical interpretation in terms of the cusp 
(eikonal) anomalous dimension; it is universal and appears in resummation of 
threshold corrections to any hard QCD process. 

However, the resummed cross section in the LLA now appears_to be sen- 
sitive to the infrared (IR) region at the level of l/Q corrections D'B. Indeed, 
remove gluons with energy Q{\ — z)/2 less than ^ ~ ^qcd by inserting the 
appropriate ^-function in the ^-integral. A simple calculation shows that the 
cross section changes then by ~ jiN/Q. Given this sensitivity to the IR region, 
one would suspect that geniune power corrections of this magnitude exist. The 
same conclusion can be obtained from considering divergences of the pertur- 
bative expansion of Eq. |^ in large orders (renormalons). One easily checks, 
however, that the dangerous IR contributions correspond to terms with less 
logarithms of N than are resummed by the LLA and thus are beyond the ac- 
curacy to which Eq. || has been derived. Thus this evidence is by itself not 
conclusive. A more accurate analysis will clarify two questions: 

• Does the IR sensitivity (of order 1/Q) of the LLA resummed cross sec- 
tion represent the 'true' magnitude of power corrections to the Drell-Yan 
process or is it artificially introduced by resummation, that is by the pro- 
cedure that separates those regions of higher order Feynman diagrams 
which give rise to large logarithms? 

• If the exact Drell-Yan cross section has no 1/Q IR contributions, can 
the resummation of large perturbative corrections to all orders be made 
consistent with the IR behaviour of finite-order Feynman diagrams? Or 
is the (spurious) l/Q sensitivity found in LLA intrinsic and unavoidable 
for resummation of threshold corrections? 

2 IR sensitivity of LLA and soft gluon emission at large angles 

We provide evidence Bill that the Drell-Yan cross section is free from 1 /Q 
IR contributions. To one-loop accuracy this statement can be checked by 
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an explicit calculation with a small gluon mass A as regulator. Nonaiialytic 
terms in the small-A expaiision correspond to higher twist contributionsQu. A 
textbook calculation givesld at 3> 1 



conviN, Q\ X^) conviN, Q\ 0) = In(AVQ^) + » ^) 

(3) 

Note that the suspected linear term -^/A^/Q^ is absent; the leading IR con- 
tribution is of order N^X^/Q^. The N'^ enhancement signals that the power 
corrections are determined by the smaller of the two large scales: Q/N is the 
moment space analogue of the energy available for gluon emission. In terms 
of the energy fraction z the relevant scale is (1 — z)Q rather the mass of the 
Drell-Yan pair Q. 

To understand the apparent 1/Q sensitivity of the LLA, it is instructive 
to consider the structure of the one-loop integral for soft gluon emission 

//l^k 
^^S[{pi+p2-kf-Q^]\Mvy\\ (4) 

The matrix element (in LLA) is |A1-dj;|^ ~ 2Q^ /k\ and it is convenient to 
rewrite the phase-space integral as 

(Pk f dkl 



2fco J 

We now take a massless gluon, and to avoid coUinear divergences introduce an 
explicit cutoff on the minimal transverse momentum k± > fi. Remembering 
that ko = \/s{l — z)/2 and taking moments we get 

The crucial point is now that the LLA corresponds to resummation of soft 
and collinear emission, that is the leading logarithms (and in fact the next- 
to-leading as well) come from the integration region of small gluon transverse 
momentum compared to its energy k± ^ ko = Q{1 — z)/2 <^ Q. Thus, for 
summation of large logarithms, it is safe to neglect the term Ak^/Q"^ under 
the square root, so that 
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where we have replaced by z^^^ — 1 to take into account virtual gluon 

exchange. Taking the integrals, we get the expected double logarithm but also 
the linear term in ^/Q discussed in Sect. 1. 

However, this IR contribution of order 1/Q comes from the end-point 
integration region 1 — z ~ ^/Q (where k± ~ kg) in which neglect of the k^/Q^ 
term under the square root is not justified. In fact, the square root cannot 
even be expanded in k'j_/{Q'^{l — z)^) since this would generate increasingly 
singulaiL contributions. Instead, the integral must be taken exactly. When this 
is doneEl, all fi/Q terms disappear. 

The physical reason for the enhanced IR sensitivity of the resummed cross 
section in LLA is that we neglected soft gluon radiation at large angles k± ^ ko. 
To recover the correct IR behavior, the phase space integral for soft gluon emis- 
sion has to be taken exactly; the common coUinear approximation is sufficient 
for summing logarithms to leading and next-to-leading logarithmic accuracy 
but is misleading for the analysis of power-suppressed effects. 

3 Resummation of soft emission and Wilson lines 

Exponentiation of large logarithms occurs for both, coUinear emission and soft 
gluon emission, including emission at large angles. Thus, after a complete 
treatment of resummation of soft gluon emission, the apparent 1/Q sensitivity 
in LLA should be compensated by other conributions to the exponent. The 
best- known generalization of the LLA formula Eq. |^ was given by Sterman eI 

\nuJDYiN,as) = / dz — <2 -rT^Tcnsp{as{k^)) 

Jo I JQHi-^-) 

+B{as{Vl~zQ)) + C(a,((l - z)Q)| + 0{l) (7) 

This expression involves three "anomalous dimensions" Fcusp,^ and C. The 
LLA corresponds to taking into account the leading term 0{as) in the expan- 
sion of Fcusp and neglecting B, C. The next-to-leading logarithmic accuracy 
requires two terms in Fcusp and the leading 0{as) terms in B, C. In general, 
higher order corrections to Fcusp j B, C generate contributions with less and less 
powers of In N for a given power of as . 

The three terms in Eq. ^ have the following origin: The function "B" 
comes from subtracting the DIS cross section (squared) and is irrelevant for 
our discussion. The term with a double integral resums soft and coUinear 
gluon radiation to all orders. This contribution is universal for all hard pro- 
cesses. Finally, "C" corrects for soft gluon radiation at large angles and is 
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process-dependent. This term starts with 0{a^) in agreement with the com- 
mon wisdom that the radiation at large angles is suppressed by two powers 
of IniV, but taking it into account can be crucial to recover the correct IR 
behavior. 

Thus, the 1/Q IR sensitivity of the generalized LLA expression given by 
the first term in Eq. |^ should be cancelled by the "C"-term. To test how 
this happens, one needs some approximation to calculate the anomalous di- 
mensions to all orders in perturbation theory. A convenient formal parameter 
is Nf, the number of light fermion flavors. The leading contribution in the 
large- A^/ limit corresponds to a chain of fermion loops inserted into the single 
gluon line. Taking into account the chain of loops has two effects: First, it 
generates the correct argument of the running coupling in Eq. ^ which is the 
gluon transverse momentum. [As usual, we tacitly assume that Nf can be used 
to reconstruct the full one- loop running, determined by Pq.] Second, counter- 
terms for individual fermion loops produce non-trivial anomalous dimensions 
Ecusp, -B, C to all orders in as- The corresponding calculation is technical and 
can be found in Ref. 1. The result is that all 1/Q IR contributions generated 
by the generalized LLA term are cancelled by the IR contributions related 
to factorial divergence of the perturbation expansion of "C". This tells us 
that although the resummation formula Eq. |7| is valid, it involves significant 
cancellations between different contributions and the true IR behaviour is re- 
stored only after summation of an infinite number of terms in the expansion of 
"C" . This contradicts the logic of resummation to resum an infinite number of 
large logarithms by calculating only a finite number of terms in the anomalous 
dimensions. 

The simplest remedy would be to use the exact phase space factor for the 
one-loop gluon emission and replace the first term in Eq. ^ by 

N-i 11 / '''^'-''>' dkl rcusp(as(fc±)) 



Jq 



(8) 

With this substitution the 1/Q IR sensitivity disappears and the functions 
B, C are modified starting 0{al) so that the undesired behavior of C in large 
orders is removed. 

A different approach to soft gluon resummation emphasizes the renormal- 
ization of Wilson lines □. Its theoretical advantage is the operator language 
that avoids the separation of small-angle and large-angle soft emission. 1/Q 
IR contributions never appear. The starting point is the well-known fact that 
soft gluon emission from a fast quark can be described by a Wilson line oper- 
ator along the classical trajectory of the quark. The product of Wilson lines 
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for the annihilating quark and antiquark is denoted by Udy (x) , where x is the 
annihilation space-time point. Up to corrections which vanish as z ^ 1, the 
partonic Drell-Yan cross section is given byO 



(9) 



^^DY = 1 + 0{as) is a short-distance dominated function, independent of z. 
Wdy is the square of the matrix element (n|rj7DY(0)|0), summed over all final 
states: 



WBY{z,as) 



dyo 

2lT 



,iyoQ{l-z)/2 ^ 



fUly{y)TUuYmO) (10) 



The Fourier transform is taken with respect to the energy of soft partons and 
y = (2/0,0). 

The crucial observation is that the Wilson line depends only on the ratio 
{liN)/{QNo) (taking moments of Wdy(z, cts)), where /x is a cutoff separating 
soft and hard emission (the renormalization scale for the Wilson line) and Nq 
is a suitable constant. Hence the iV-dependence of the Drell-Yan cross section 
in the soft limit can be obtained from the /i-dependence of Wilson lines, which 
is given by the renormalization group equation (here Us = c^si^J■)) 



2A 



_d_ 

da. 



InV^DY 



, a. 



= r. 



cusp 



(as) In 



DY 



(11) 

It involves two anomalous dimensions rcusp(as) and rDY(as) related to the 
cusp and to presence of light-like segments on the Wilson line, respectively. 
The general solution of Eq. |ll| is given by 



IuIVdy 



as{Q) =lnT4^DY(l,a.(QiVo/iV))+ 



rcusp(as(fc±) In 



+ TuYiasik^)) 



(12) 



where we set fi — Q. The inhomogeneous second term in Eq. |lj can be 
rewritten (identically) in a more familiar form, which resembles the first term 
in Eq. 0: 



Jo 1-2 Jr 



dk"^ 

— j"rcusp(as(fc±)) + TY)Y{as{{l - z)Q) 



(13) 
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Note presence of the initial condition W-Dvi^, cts{QNQ/N)). Its expansion in 
as produces subdominant logarithms a^ln''^^ A'^ which can be absorbed into 
a redefinition of Fdy ■ 

rDY(«.) -^C(as) =rDY(as)-/?(a,)-^ lnWDY(l,a.). (14) 

das 

C{as) starts at order a^. It does not affect resumniation of large logarithms 
in N to next-to-leading accuracy a^ In^ N. 

It remains to subtract the DIS cross section, which can also be implemented 
in the language of Wilson lines, see Ref. 9 for details. Finally, we get 

rl-No/N ^ ( rQ^il-zf 

lncjDY(A^,as) = -/ dz- {2 ^r,^,p{asik^)) 

Jo I "'q^(i-z) ^-L 



+ B{as{VT^Q)) + C{as{{l - z)Qj) } + 0(1) . (15) 



(A^o = exp(— 7£;) in the MS scheme.) This form of the resummed cross section 
is as legitimate in the framework of the perturbation theory as the more con- 
ventional expression in Eq. |^. They have different properties, however, as far 
as sensitivity to the IR behavior of the coupling is concerned, which becomes 
important when the anomalous dimensions Fcusp, ■ • ■ are truncated to finite or- 



der. Since the region of very large z ^ 1 is removed in Eq. 15, this expression 
shows no IR sensitivity at all unless N > Q/Aqcd- Loosely speaking, this is so 
because the Wilson line approach treats small- and large-angle gluon emission 
in a coherent way. Because this technique can also treat subleading logarithms 
in a systematic way, it is preferred over, for example, the modification of the 
phase space as in Eq. |^. p. 

It is natural to expect (and explicit calculationld in the large- iV/ limit con- 
firms this) that the anomalous dimensions Fcusp (c^s) and Fdy(q!s) in the MS 
scheme are analytic functions of the coupling at as — 0. Then, all power cor- 
rections to the resummed cross section (to all orders in NAqcb/Q) originate 
exclusively from the initial condition for the evolution equation for the Wilson 
line, and are not created (or modified) by the evolution, i.e. by soft gluon re- 
summation. Thus, if the resummation of soft gluon emission is done coherently 
for all angles, the only effect of soft gluons on power corrections is a change 
of scale, the replacement Q — > Q/N as the parameter of the power expansion. 
This suggests that, in general, there is no reason to suspect new nonpertur- 
bative contributions in resummed cross sections as compared to finite-order 
calculations. The conclusion that power corrections to Drell-Yan production 
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are suppressed by 1/Q^ is then coiiaistent with the analysis of power corrections 
at tree level by Qiu and StermanE2l. 

The redefinition in Eq. |l^ transforms the IR sensitivity (and potential 
power corrections) of the initial condition for the evolution equation for the 
Wilson line into IR sensitivity of the function C in Eq. |l5| As mentioned 
above, this function becomes important precisely when one starts to be sensi- 
tive to gluon radiation at large angles, and the conclusions on power correc- 
tions depend sensitively on this region. Because of this, we are sceptical that 
universality of nonperturbative corrections to resummed cross sections could 
be deduced from the universality of soft-coUinear gluon emission as embodied 
by the eikonal (cusp) anomalous dimension, an idea originally put forward in 
Refs. 3,11. In the Wilson line technique the solution of the evolution equation 
never involves the QCD coupling integrated over the IR Landau pole as long 
as < Q/Aqcd- Since this inequality sets the boundary for a perturbative 
treatment anyway, Eq. |l^ (which coincides with the resummation procedure 
used in Ref. 12 to next-to-leading logarithmic accuracy) is suited for all mo- 
ments that can be treated in a power expansion in 1/Q. 

4 Top quark production at the TEVATRON 

In the light of our discussion, let us consider recent results for the resummed top 
quark production cross section, which we sumraajjize in Table 1. We concen- 
trate on the comparison of two new calculationsEjO. Both assumed = 175 
GeV and used the same parametrisations for the structure functions. Hence 
the difference is entirely due to different resummation prescriptions. The dif- 
ference in central values is of order 15%, compared to ^ 10% renormalization 
scale dependence and 5% due to uncertainty in the structure functions. Ap- 
parently, resummation causes the largest ambiguity. Note that the resummed 
cross section of Ref. 14 practically coincides with the strict 0{a^) result. Thus, 
in Ref. 14 the resummation of In A^ terms has a negligible effect, while the re- 
summation in Ref. 13 produces a 10—15% enhancement. Since both procedures 
sum all leading logarithms (in the sense of Eq. |^) , the difference is entirely due 
to terms with less powers of logarithms which are beyond the accuracy of the 
resummation in the strict sence. Unless we can prefer one particular resum- 
mation procedure, the difference would have to be considered as the present 
theoretical uncertainty. Our discussion of Drell-Yan production suggests the 
criterium that resummation procedures should not introduce power correc- 
tions (factorial divergence in large orders) which are not already present in 
finite order approximations. From this point of view, we are led to prefer the 
prescription used in Ref. 14, which starts from Eq. Oa. 
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Table 1: Resummed cross section for the tt production at the TEVATRON. mt = 175 GeV, 
MRSA' parton distributions for the central value quoted. 



Rcf. 




Uncertainty 


BC 13 
CMNT 


\3 


4.95 
5.32 
4.75 


4.53 - 5.66 
4.93 - 5.40 
4.25 - 5.00 



The major numerical difference between Ref. 13 and Ref. 14 comes from 
a different procedure to implement the inverse Mellin transformation from 
moment space to momentum space. The subtle problem here is to which ex- 
tent one can avoid contributions of very large moments N > Q/Aqcd, which 
strictly speaking can not be treated by short-distance methods. This problem 
is somewhat similar to the problem of analytically continuing perturbative 
QCD predictions from Euclidian in Minkowski space, relevant, for example, in 
connection with the r-lepton hadronic width. The particular way of perform- 
ing the analytic continuation becomes important when one uses a resummed 
coupling constant, and the guiding principle proves to be to avoid the region 
Q < Aqcd in the complex Q-plane, where no short-distance treatment is pos- 
sible. If the region Q < Aqcd is not avoided, one may introduce spurious l/Q^ 
power corrections tathe decay widthQ, which are absent in the OPE. Similarly, 
Catani et al. find E-3 that the inverse Mellin transform of the resummed cross 
section in moment space has to be done by exact numerical integration in the 
complex N plane, avoiding the region Re N oo where the IR singularity 
in the running coupling becomes important. Failure to avoid this region may 
result in spurious effects of order {Aqcd/Q)^"'^ ■ 

Within the approach of Ref. 13 this problem is somewhat masked by us- 
ing a resummation formula similar to Eq. 0, in which the sensitivity to the 
IR behavior of the coupling is of order 1/Q for any N. As explained above, 
this IR sensitivity is an artifact of truncating the anomalous dimensions to 
finite order. One also notes that applying the prescription of Ref. 13 to DrcU- 
Yan production requires to introduce a phenomenological 1 / Q-correction, that 
mainh^ seems to cancel the l/Q-effects generated by the resummation prescrip- 
tion lil. 

The difference between various resummation procedures should also be 
perceptible in high-p^ jet production at the TEVATRON, and a theoretical 
understanding of this difference might be one aspect in understanding the 
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apparent excess of large-pj_ jets seen by CDF. 
5 1/Q IR sensitivity of thrust 

The Drell-Yan cross section appears to have no l/Q power corrections. This is 
not generally the case for any quantity. There are good reasons to suspect the 
existence of l/Q nonperturbative_hadronization corrections to event shapes 
observables in e+e~ annihilation EZl. Unlike the Drell-Yan cross section, these 
observables cannot be expressed directly in terms of Feynman diagrams, and 
are obtained by integrating the QCD amplitudes with certain weight functions 
such as to emphasize a particular final state configuration. These weight func- 
tions generally destroy the balance of gluon emission at small and large angles, 
and make these observables sensitive to nonperturbative momentum flow at 
large angles. As a consequence l/Q corrections are invariably expected for 
event shapes. For example, the average thrust (1 — T) of the final state is 
computed to leading order in by inserting the factor 

1 - T = (fco - ^kl-kl)/Q (16) 

into the phase-space integral for gluon emission, which has a structure identical 
to the Drell-Yan cross section in moment space, see (^. The above factor 
suppresses small-angle emission but causes l/Q IR sensititvity. 

An interesting speculation is whether 1/Q corrections to event shapes are 
universal ip-the sense that they can be related to a single nonperturbative 
parameter OJOEJ. Although, due to importance of large angle emission, this 
parameter would not be related to the universal cusp anomalous dimension, 
the hypothesis makes sense so long as the underlying physical process is the 
same for all event shapes. Strictly speaking, the answer sepps to be negative, 
since 1/Q corrections also occur outside the two-jet regiorc3 and higher-order 
corrections are not suppressed, because is evaluated at low scale, so thaLit 
is not counted. One may still argue in favour of approximate universalityQci, 
if the coupling stays finite and reasonably small in the infrared. This purely 
phenomenological hypothesis could in principle be subjected to experimental 
tests. This, in fact, seems very hard in practice, because of poor control of 
higher orders of the perturbative series. One may suspect that the largest 
part of the hadronization correction to event shapes estimated by Monte Carlo 
event generators is in fact related to the perturbative parton cascade which can 
indeed be universal to the extent that the event shape variable is dominated 
by the two-jet kinematics. 

To illustrate the difficulty in testing universality consider the average 
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thrust (1 — T). The existing experimental data are weU described by0 

(l-r>(Q) = 0.335 a,(Q) + 1.02 a,(g)2 + i^ (17) 

where the first two terms give the QCD calculation (to O(a^) accuracy). The 
power correction is needed to gain agreement with the data (over the entire 
range of Q). The second order perturbative correction has a large coefficient, 
indicating that the adopted scale Q is in fact inadequate for this process. The 
scale setting problem for event shapes is difficult. However, as a natural first 
guess one can try to take as at a scale of order of the jet mass M, which 
is related to thrust in the two-jet fimit by = (1 - T) Q^/2 [We take the 
scale Ql = {1 — T) Q^/A, since for fixed T this is the upper limit on the gluon 
transverse momentum.]. Taking as{mz) = 0.12 and using the fixed-order QCD 
result (1 — T) ~ 0.07, we get ~ O.ISQ. Fitting again a l/Q correction to 
the same data, we get 

(l~T)(g) = 0.335a,(g,) + 0.19a.(Q*)' + (18) 

The second-order coefficient in the perturbative series has become much smaller 
and the size of the required hadronization correction has also decreased. One 
might actually think of writing the power correction as O.OSGeV^/Q*. Viewed 
this way, the issue of universality becomes inseparable from the problem of 
determining the most appropriate scale for the process. 
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